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SUMMARY 

A viscoplastic constitutive theory is presented for representing the high- 
temperature deformation behavior of metal matrix composites. The point of view 
taken is one of a continuum wherein the composite is considered a material in 
its own right, with its own properties that can be determined for the composite 
as a whole. It Is presumed that a single preferential (fiber) direction is 
identifiable at each material point (continuum element), thereby admitting the 
idealization of local transverse isotropy. A key ingredient in this work is 
the specification of an experimental program for the complete determination of 
the material functions and parameters for characterizing a particular metal 
matrix composite. The parameters relating to the strength of anistropy can be 
determined through tension and torsion tests on longitudinally and circumferen- 
tially reinforced thin-walled tubes. Fundamental aspects of the theory are 
explored through a geometric interpretation of some basic features analogous 
to those of the classical theory of plasticity. 


INTRODUCTION 

Structural alloys used in high-temperature applications exhibit complex 
thermomechanical behavior that is Inherently time-dependent and hereditary, in 
the sense that current behavior depends not only on current conditions but also 
on thermomechanical history. Considerable attention is being focused now on 
metal matrix composite materials that possess strong directional characteris- 
tics. In high-temperature applications these materials exhibit all the com- 
plexities of conventional alloys (e.g., creep, relaxation, recovery, rate 
sensitivity, etc.), and in addition, their strong initial anistropy adds fur- 
ther complexities. 

Here, we present a continuum theory to represent the high-temperature, 
time-dependent, hereditary deformation behavior of materials that are initially 
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transversely Isotropic. The theory is intended to apply to materials, particu- 
larly metallic composites, that can be idealized as psuedohomogeneous continua 
with locally definable directional characteristics. 

The composite material is viewed as a material in its own right, with its 
own properties that can be measured and specified for the composite as a whole. 
Experiments for this purpose are outlined in detail in the fourth section of 
the paper. This view is intended to satisfy the structural analyst or design 
engineer who needs reasonably simple continuum methods of structural analysis 
to predict deformation behavior in complex multiaxial situations, particularly 
at high temperature where material response is enormously complex. Indeed, the 
prediction of component lifetime depends critically on the accurate prediction 
of deformation behavior. 

The alternative approach is concerned with detailed Interactions of the 
constituents of the composite: fabrication, bonding, and the relation of the 

properties of the composite to the Individual properties of the fiber and 
matrix. Clearly, such problems are of great importance, and the two approaches 
mentioned are not mutually exclusive. Here, however, the continuum point of 
view will be emphasized. This is done in the same spirit that the theories of 
elasticity, plasticity, viscoelasticity, and others are formulated; on the 
basis of macroscopic observations, without direct consideration of the details 
of i ntermolecular , intergranular or interdislocation Interactions. Of course, 
this is not to imply that qualitative (and quantitative) understanding of 
behavior on the microscale should not strongly influence the formulation and 
structure of phenomenological theories. 

The authors are hopeful this research will complement other ongoing 
efforts at NASA Lewis Research Center relating to the high-temperature behavior 
of metal matrix composites (refs. 1 to 5). Parts of the present work (believed 
essential in representing the time-dependent, hereditary behavior of metals) 
may prove helpful in extending the micromechanics equations for the thermal 
and mechanical behavior of composites (refs. 1 and 2) to Include some important 
viscoplastic features. 


STATEMENT OF THE THEORY 

This work is an extension of that by Robinson (ref. 6) and includes the 
former work as a special case. In reference 6, three material parameters over 
and above those necessary for representing isotropic viscoplastic behavior were 
necessary to account for transverse isotropy. Here, four parameters have to be 
specified (fig. 1). The additional parameter arises from a less restrictive 
set of assumptions made in the theoretical development. Definition of the 
additional parameter leads to more testing to characterize a particular material 
but, at the same time, offers the distinct advantage of greater flexibility in 
correlating predictions with experimental data. 

As in reference 6, the starting point is the assumed existence of a dis- 
sipation potential function (refs. 7 to 9); that is 

Q =» Q (a. j ,a. j ,T) (1 ) 
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with the generalized normality structure 
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Here, a.^ and a^ denote the components of the applied and internal stress 
tensors .respectively, denotes the components of the inelastic strain rate 
tensor, h is a scalar function of the internal stress, and T is the tempera- 
ture. Although Q is shown as a function of temperature, only isothermal 
deformations will be considered in the following development. Extension to 
noni sothermal conditions follows as in reference 10 . 


In the fully isotropic case, the stress dependence of Q enters only 
through the principal invariants of the deviatoric applied and internal 
stresses (ref. 10). For transverse isotropy, 0 must depend additionally on 
the local preferential (fiber) direction denoted by the components of a unit 
vector d, (or, as the sense of d, is immaterial, on the components of a 
symmetric directional tensor d^dj). Form invariance (objectivity) of Q 
requires that it depend only on invariants of the applied and internal 
stresses, the directional tensor, and certain products of these tensors 
(ref. 11 ). 


A subset of the irreducible set of invariants for form invariance (integ- 
rity basis) is used (ref. 11 ) 
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where denotes the components of the effective stress, that is, the dif- 
ference of the deviatoric applied and internal stresses. I] relates, as in 
the isotropic case, to the effective octahedral shear stress, I 2 relates to 
the shear component of the effective traction on the plane of isotropy (plane 
normal to d^), and I 3 corresponds to the normal component of the effective 
traction along dj . 

Taking Q to be dependent on the appropriate Invariants and using equa- 
tions (2) and (3), the flow law becomes 
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and the evolutionary law becomes 
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and 
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The function f(F) and the material parameters m, 6 , R, H, and Ky are asso- 
ciated with the viscoplastic response (ref. 10 ); the components of the unit 
vector di (specified in terms of two Euler angles), u, and n. are the four 
parameters associated with the direction and strength of transverse isotropy 
(fig. 1 ). I ^ , I ^ , I 3 are invariants of the deviatoric internal stress a..., 

similar in form to equations (4) to (6). Note that with co = n = 1 
equations (7) to (10) reduce to the isotropic forms reported in reference 10. 
The details of the derivation of equations (7) to (10) are left to the 
references 6 , 10 , and 1 1 . 


For a particular composite material the parameters u and n» designat- 
ing the strength of anisotropy, depend on the individual constituent materials 
(fiber and matrix) and their volume ratio. Different volume ratios involving 
the same constituents are considered different materials. In extending the 
present theory to arbitrarily large deformations, the local volume ratio may 
change in the course of deformation (as does the local density in an isotropic 
material); in which case evolutionary laws for w and n must be specified. 
Also, for large deformations and/or rotations, the preferential direction d^ 
may convect with the material thereby resulting in increasing anisotropic inho- 
mogeneity. At this time, large deformations and rotations are not considered. 
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SOME FUNDAMENTAL IMPLICATIONS OF THE THEORY 


As in earlier works (refs. 6 and 10) F plays the role of a Blngham-Prager 
threshold stress function; inelastic response occurs only for F > 0. The sur- 
face F = 0 in the stress space encloses stress states that produce elastic 
behavior only. Figure 2 shows a typical projection of F = 0 on the - o 22 

plane for full isotropy (u * n = 1) and the virgin state (a|j = 0). An infi- 
nite family of surfaces F = constant is associated with eacn Inelastic state. 

The direction of the inelastic strain rate vector at each stress point on a 
given surface is directed normal to the surface. The existence of these sur- 
faces and the concept of normality has been demonstrated experimentally for the 
Isotropic case in reference 12. 

Figure 3 shows the corresponding projection of F = 0 on the - a 22 

plane for the transversely isotropic case with w > r» * 2. For each curve 
shown, the preferential direction is taken to lie in the X] , X 2 physical 
plane with a specified angle <p relative to the X] axis. These curves, as 
in figure 2, correspond to the virgin state. The shape and orientation of the 
surfaces F = 0 (and all surfaces F - constant) now depend on the local pref- 
erential direction di . Note, for Instance, that for <p - 0, the Intercept on 
the o-j j axis Is the threshold stress Y|_, (figs. 1 and 3), the Intercept on 

the c 22 axis Is Yj, and w = (Y(_/Yj) - 2. 

As an Interesting and Illustrative example, consider the stress path 
°11 “ 2o 22 1s denoted as a dotted line In figures 2 ahd 3. This Is 

equivalent to the stress state In a thin-walled tube with closed ends under 
Internal pressure, where < 7 ^ Is the circumferential or hoop stress and o 22 

the axial stress. With increasing pressure, the stress point eventually 
reaches the surface F = 0 and the tube begins to deform 1 ne 1 as 1 1 cal ly. In fig- 
ure 2, corresponding to the Isotropic tube, Inelastic deformation occurs as 
Indicated by the strain rate vector shown. Normality dictates that the axial 
strain rate e 22 is zero, that Is, the tube Incurs no Inelastic change In 
length. Contrast this behavior with that of figure 3 with <p * 0. This case 
represents a circumferentially reinforced tube with a threshold stress In the 
circumferential direction that is twice that of the axial direction. As F » 

0 Is reached and inelasticity begins to <pccur, the (normal) strain rate vector 
has a relatively large axial component e 22 . The thin tube now experiences 

inelastic axial extension. Thus the mode of inelastic deformation has changed 
qualitatively with reinforcement. Similar observations are well documented 
for time-independent reinforced structures (ref. 13). 


EXPERIMENTAL DETERMINATION OF MATERIAL PARAMETERS 

Two types of specimens are presumed to be available: thin-walled compos- 

ite tubes that are longitudinally reinforced (having a single fiber direction 
oriented axially) and those that are circumferentially reinforced (circumferen- 
tial fiber orientation). Each type of tube will be loaded either In pure tor- 
sion or in pure tension. Although not discussed here, combined tension and 
torsion experiments can be used as verification tests to assess the correct- 
ness of the multiaxial theory (ref. 12). As is well known, the thin-walled 
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tube Is an ideal specimen for the development of constitutive relationships in 
that it provides a nearly homogeneous and uniform region of stress and strain, 
and is statically determinate. 

Those parameters relating to the strength of anisotropy, w and n, and 
the threshold strength Y|_(=Kj(4«2 - 1)1'2) are determined through probing 
tests. These tests are designed to determine the inelastic strain rate for a 
given stress in the neighborhood of a constant inelastic state, here the ini- 
tial (virgin) state of the material. Indeed, it is the degree of initial 
transverse isotropy that is of interest in this study. The probing test, con- 
ducted properly, furnishes the desired information without significantly chang- 
ing the state. 

The material function f(F) and the parameters m, (3, R, and H relat- 
ing to the viscoplastic properties of the composite are obtained from a combi- 
nation of the probing tests and creep tests conducted on longitudinally- 
reinforced tubes (or bars). Alternately, in place of the latter tests, one 
could use torsion tests on circumferentially reinforced tubes as outlined in 
reference 11. The present choice of basing the viscoplastic parameters on uni- 
axial creep tests of axially reinforced tubes is motivated by (1) the relative 
ease of fabricating longitudinally reinforced tubular specimens over those 
reinforced circumferentially and (2) the advantage of characterizing directly 
the inelastic response in the critical fiber direction. In the case of extreme 
reinforcement (e.g., a relatively high volume ratio of very strong fibers that 
remain essentially elastic), composite structures are known to be "shear limit- 
ed" (ref. 12), and their inelastic behavior is governed largely by the shear 
response of the matrix. Under these conditions it may be advantageous to 
determine the viscoplastic parameters through the torsional creep tests on cir- 
cumferentially reinforced tubes discussed in reference 11. 

First, consider a probing test on a longitudinal tube. Pure torsion and 
pure tension probes are schematically Illustrated as the respective paths o-a 
and o-b in the a-x stress space of figure 4(a). Data from such tests take 
the form of a sequence of stress and inelastic strain rate pairs, (x,y) along 
o-a and (o,c) along o-b. These data can be conveniently plotted as the solid 
curves o versus e and x versus y in figure 5. Extrapolation of the 
o versus e curve in figure 5 to the e = 0 axis furnishes the longitudinal 
threshold stress Y[_. 

Specialization of equations (9) and (7) for the path o-a in figure 4(a) 
results in 


F = 


and 


in which x(=a 12 ) is the applied shear stress. The corresponding equations 
for path o-b are 
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where cKEa^) Is the applied normal stress. Now for F = constant (Q = con- 
stant), equations (11) and (13) give 



while equations (12) and (14) combine to give 



(15) 
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or 


oe » zy ( 1 7 ) 

Thus, In figure 4(a) points (t,y> along o-a and (o,e) along o-b, having 
the same dissipation rate, lie on a common F = constant (Q = constant) curve; 
for example, the particular points (-c A> y A > and (o A ,e A ) In figure 4(b). 

Pairs of points In the plot of figure (5) that lie on an F = constant 
locus are related geometrically such that areas and t a y a are equal. 

Several such pairs can be matched up giving an average value of the ratio 



Thus, from equation (15) 



(18) 


(19) 


Probing tests on circumferentially reinforced tubes produce results 
entirely analogous to those discussed above for longitudinally reinforced 
tubes; counterparts of figures 4(a), 4(b), and 5 can be constructed. 


The governing equations corresponding to the pure torsional loading of 
the circumferential tube are Identical to equations (11) and (12). The equa- 
tions relating to pure tension of the circumferential tube are 
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c n = e = f(F) —I 2 — to 2 (21) 
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As before, for F = constant (Q = constant), equations (11) and (20) combine to 
give 
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whereas equations (12) and (21) give 
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Again, In figure 4(a) points (x,y) along o-a and (o.e) along o-b, with equal 
dissipation rates, fall on a common F = constant curve (fig. 4(b)). Simi- 
larly, matching pairs of points corresponding to F - constant curves and 
averaging gives 



(25) 


Making use of equation (19), from tests on longitudinally reinforced tubes, and 
equation (25), from circumferentially reinforced tubes, results In 


and 
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thereby completely characterizing the strength of initial anisotropy. Recall 
that the longitudinal threshold stress Y|_(=Kj(4u2 - 1 )*' 2 ) a i so known. 

Turn now to the determination of the remaining viscoplastic material 
parameters. As indicated earlier, although several options are open in this 
regard (including torsional tests on circumferentially strengthened tubes as. 
in ref. 11), the choice here is to consider uniaxial creep tests on longitudi- 
nally reinforced thin-walled tubes. This is in addition to the probing tests 
on longitudinal tubes already discussed. Typical results of uniaxial creep 
tests are illustrated in figure 6. 

According to the present theory, the governing equations for the consid- 
ered creep conditions are 


and 


where 
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Here, etEe^) is the axial component of inelastic strain rate, and s(=a 1 |) Is 
the uniaxial component of the internal state variable . The first of equa- 
tions (32) indicates that the function f(F) has been specialized as a power 
function characterized by the constants p and n. 
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Note for future reference that for F >> 0, equation (30), the flow law, 
can be approximated as 


c = B(a - s) 


2n+l 


where 


y(4u 2 - 1 )Y 2n 


(33) 


(34) 


Also, during the early part of primary (transient) creep where s is small 
(fig. 6), the evolutionary law equation (31) can be approximated as 


s 



(35) 


The initial creep rate following abrupt application of the stress o 
(fig. 6), is expressed by equation (33) with s * 0. That is 
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Bo' 


2n+l 


(36) 


Information for determining B and n can be obtained directly from the ini- 
tial creep rates as Illustrated in figure 6, or alternatively, and more accu- 
rately, from the data generated from the probing tests already considered; the 
o versus e data Illustrated in figure 5. By assuming these data correlate 
with equation (36), a plot of log(e) versus log(o) provides both n and B 
directly (fig. 7). If correlation with equation (36) Is not satisfactory a 
different function f(F) may have to be considered in the first of equa- 
tions (32). In this study, it is assumed that the power law form is appropri- 
ate (as has been found for several isotropic alloys) and that n and B can 
be determined. With n, B, to, and Y[_ known, y is then determined from equa- 
tion (34). 


Now focus attention on steady state creep information (fig. 6). At 
steady state s = 0, so from equation (31) 


c 


s 



(37) 


Steady state creep data provides the pairs (a,e s ) but with w, y, Y|_, and n 
known, equation (30) allows s s (the steady state internal stress) to be calcu- 
lated for each pair, thus giving the data pairs (e s ,s s ). Plotting log(e s ) ver- 

/N /N 

sus log(s s ) provides values of ( R/H ) and m directly as Indicated by the 
logarithmic form of equation (37) and figure 8. The last of equations (32) 
provides the exponent m. 
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Attention is now turned to the primary creep stage of the creep tests. As 
indicated in equation (35), the evolution of s in the early stages of primary 
creep is governed by 


(38) 



Eliminating time and Integrating 
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s 13 ds 



results in 


(39) 
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or in transposed logarithmic form 
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Now using equation (33) to solve for s yields, 
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a 



l/2n+l 


(42) 


and substituting this result into equation (41) gives 

+ log ) (43) 

VH(3 + 1)/ 

Early primary creep data provides the data triplets (a,e,e) at each time. With 
B and n known, these data can be plotted (fig. 9) in the form of equa- 
tion (43) yielding H(B + 1) and (6 + 1). Using equations (32), 6 and H 
can then be determined. Since R/H is known from steady state creep data, 

/N 

individual values for R and H can then be found. Further, by making use 
of the second and third of equations (32), the parameters R and H then are 
known. Finally, the viscoplastic parameters Yl(EKj(4u 2 _ l) 1 ' 2 ), ^ n, ^ 
R, and H are completely determined, as are the measures of the strength of 
transverse isotropy, to and n- 


log(e) = (3 + 1 ) log 



l/2n+l, 


SUMMARY AND CONCLUSIONS 

A constitutive theory is presented to represent the high-temperature, 
time-dependent, hereditary behavior of materials that can be idealized as 
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initially transversely Isotropic. In particular, the theory is applicable to 
metal matrix composite materials at elevated temperature where their mechanical 
behavior includes significant viscoplasticity (e.g., creep, relaxation, thermal 
recovery, etc.) and, at the same time these materials are strongly directional. 
It is presumed that a single preferential (fiber) direction is identifiable at 
each material point thereby admitting the idealization of local transverse iso- 
tropy. Although not addressed here, the theory can be extended, at the expense 
of some additional complexity, to account for two (or more) identifiable pref- 
erential directions at each material point. 

The composite is viewed as a continuum in its own right; and detailed 
interactive effects of the constituents are not accounted for directly. Of 
course, this precludes predictions of detailed phenomena such as failure by 
debonding, delamination, and so forth. However, the result is a reasonably 
simple multiaxial constitutive theory that is easily implemented into struc- 
tural analysis codes for predicting the deformation response of structures sub- 
jected to complex thermomechanical loading histories. Because the response in 
the presence of material anisotropy is often highly nonintuiti ve, this theory 
provides a valuable tool for the design engineer. 

Some fundamental aspects of the theory are explored through geometric 
interpretation of some basic features analogous to those of time-independent 
plasticity theory. Convexity of the dissipation potential surfaces (F = 
constant or Q = constant) is demonstrated, and the shape of the surfaces is 
shown to be dependent on the strength and orientation of anisotropy. An exam- 
ple involving the response of a thin-walled tube under internal pressure demon- 
strates the qualitative changes in the inelastic deformation mode that can 
result from directional strengthening (anisotropy). 

A key ingredient in the present work is the specification of an experi- 
mental procedure for the complete determination of the material parameters for 
a particular metal matrix composite. The parameters relating to the strength 
of anisotropy are determined through probing experiments on thin-walled tubes 
of two kinds; circumferentially reinforced (a single fiber direction oriented 
circumferentially) and longitudinally reinforced (axial fiber direction). The 
tubes are loaded in both tension and torsion. The parameters relating to the 
viscoplastic properties of the composite are determined primarily through uni- 
axial creep tests conducted on longitudinally reinforced tubes. Alternately, 
as discussed, one could use pure torsion tests on circumferentially reinforced 
tubes in place of the uniaxial creep tests. Additional tests are suggested in 
order to assess the correctness and accuracy of the theory. 
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« = y l /y t n - k l /k t 


Figure 1. - Parameters accounting for transverse isotropy 0, 
u , rj . The local fiber direction is denoted by the unit vector 
dj - SPECIFIED BY THE EULER ANGLES l| / AND 0. Y^ Y y , AND 

ARE THRESHOLD STRESSES AS INDICATED* DEFINITIONS OF THE PARAM- 
ETERS U AND Q ARE GIVEN. 



Figure 2. - Projection of flow surface F = 0 on o^~o 2 2 plane 
( o = H “ 1 - i sotropy ) . Units of stress are arbitrary. 
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Figure 3. - Projection of flow surfaces F = 0 on PLANE 

(o = n = 2). Units of stress are arbitrary. 
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Figure 4(a). - Tension and torsion probes schematically illustrated 

IN THE 0-T STRESS PLANE. 

4(b) . - Data points that lie on a common F = constant curve. 
Points A, A have the same dissipation rate. 
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